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ABSTRACT 

We use the theory of Yennie, Frautschi and Suura to realize, via Monte Carlo 
methods, the process / / — > f'f' + rry at SSC and LHC energies, where / and 
/' are quarks or leptons. QED infrared divergences are canceled to all orders in 
perturbation theory. The resulting Monte Carlo event generator, SSC-YFS2, is 
used to study the effects of initial-state photon radiation on these processes in the 
SSC environment. Sample Monte Carlo data are presented and discussed. We find 
that the respective multiple-photon effects must be taken into account in discussing 
precise predictions for SSC physics processes. 
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1. Introduction 



Now that the SSC is under construction, it is extremely important to prepare 
for the maximal physics utilization and exploration of the new frontier which it 
will probe. In particular, it is of some import to determine the effects of higher- 
order radiative corrections on the SSC physics processes of interest so that optimal 
discrimination between signal and background can be realized. In this paper, we 
explore the first step in the determination of such corrections by computing them 
for the basic QED-related effects in f f -> f'f) + n-y at SSC energies using 
the methods which two of us (S.J. and B.F.L.W.) introduced for the analogous 
processes e + e~ — > // + wy for high precision Z° physics at SLC and LEP. Thus, 
here, we extend the SLC/LEP Monte Carlo event generator YFS2 Fortran in the 
first paper in Ref. [1] to the SSC physics environment. An analogous study of 
the multiple-gluon radiation in SSC processes such as qq — > q'q 1 + nG (where q, q' 
represent quarks and G a gluon) will appear elsewhere [2]. 

In the SSC environment, multiple-photon and multiple-gluon radiative effects 
are expected to be important, in partial analogy with the significance of multiple- 
photon radiation in the SLC and LEP environments. The analogy is partial because 
the would-be resonance at the SSC, the Higgs, is actually quite broad compared 
to the Z° at SLC and LEP. Nonetheless, if k « 0.001^/5/2 is a typical infrared 
resolution factor for an SSC detector, then the probability of an incoming -u-quark 
to radiate at the SSC is, e.g., 

P(ko < k < V~s/2) » ^(log(^/6m u ) 2 - 1) log(V5/*o) n n 

71 (1.1) 

= 0.44 , 

where y/s = 40 TeV and m u ~ 5.1 MeV. Hence, such radiation and its attendant 
effect on the respective SSC event structure must be computed to the standard 
SLC/LEP precision to assess its interplay with detector cuts, physics signals and 
physics backgrounds . This would then leave only the multiple-gluon radiative ef- 
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fects to be taken into account to gain a complete view of higher-order radiative ef- 
fects to SSC physics processes. Such gluon radiation will be taken up elsewhere [2]. 

Our strategy is to treat the incoming quark radiation via the YFS theory [3] 
so that we realize it on an event-by-event basis using the methods in Ref. [1]. 
The full multiple-photon character of the final state, including the physical four- 
momentum vectors of the photons, is then made available to the users of the 
attendant new multiple-photon event generator SSC-YFS2. The implementation of 
arbitrary detector cuts on the respective simulated cross sections is straightforward. 

A logical next step is to include the effects of the final-state multiple-photon 
radiation via the extension of the Monte Carlo event generators BHLUMI2.0 For- 
tran [1,4] and YFS3 Fortran [3] to the SSC environment in analogy with our ex- 
tension of YFS2 Fortran in the current work. We shall discuss these extensions in 
a future publication [2] . 

Our work is organized as follows in the paper: in the next Section, we review 
the YFS methods as they are implemented in YFS2 in the first paper in Ref. [1], so 
that this paper is self-contained; in Section 3, we discuss how we extend YFS2 to 
SSC processes and energies to get the Fortran Monte Carlo event generator SSC- 
YFS2; in Section 4, we present some sample Monte Carlo data for SSC physics 
processes and comment on their implications; finally, in Section 5, we present our 
outlook and summary remarks. 

2. Review of YFS methods 

In this Section, we review the methods used in Ref. [1] to realize the YFS 
theory via the Monte Carlo event generator YFS2 Fortran for e + e~ — > // + 717, 
/ ^ e, in the Z° energy regime. We begin by recalling the key ingredients of these 
methods. 

The YFS Monte Carlo methods in Ref. [1] take advantage of the expansion of 
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the total cross section for the process illustrated in Fig. 1, 



e + (pi) + e-(p2) — f(qi) + /fe) + 7(^1) + ■ ■ ■ + l{k n ) , (2.1) 

in terms of the YFS hard-photon residuals [3] f3 n (k\, . . . , k n ), which are free of all 
virtual and real infrared divergences to all orders in the QED coupling constant a, 
and of the products of the YFS infrared emission factors [3] 

§(*) = -« (*__*>)'+... , (2.2) 
4tt 2 \pi ■ k p 2 ■ k) 

where . . . represents the remaining terms obtained from that shown by the appro- 
priate substitutions of {(e,pi) , (— e,p 2 )} with {(Qfe,q\) , (— Qfe, (72)}, with due 
attention to signs associated with the direction of the flow of charge. The most 
infrared-singular contribution to the cross section involves n factors of S, as we see 
from the following expression for the respective cross section: 

daW = (S(h) ■ ■ ■ S(k n )P (p h p2, gi, ffi) + ... + P n (h, kn)j 

1 4 / A \ Al d 3 © d% Sk n 2aKeB ^ 



X 



where B is the YFS virtual infrared function and is given in Refs. [1] and [3]. 
Basing ourselves on Eq. (2.3), in YFS2 Fortran we proceed as follows. 

We use the YFS form factor, 



(2.4) 



(f d 3 k ~ 
2aRe£+ / -^-S(pi,p2, k) (1 - 9(k° - ey/s/2)) 

= exp (^-(2(\n(s/m 2 e ) - 1) lne + \ \n(s/m 2 e ) - 1 + n 2 /3^j , 

to compensate for the omission of small-energy photons with k° < e^fs/2 for 
e <C 1 from the phase space in (2.3) to all orders in a, as is effected by inserting 
nT=l Q^ki/y/s — e) into (2.3) and summing over all AtH Here we presume we 
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are in the e + e~ center-of-mass frame. For YFS2, the hard photon residuals /?o 1 2 
are used, as two of us have explained in Ref. [1]. This means that, in the YFS2 
Monte Carlo itself, some choice must be made for the reduction of the n-photon 
+// phase space to the j-photon +// phase space (n = 0, 1,2, . . ., j = 0, 1,2, 
n > j), which is involved in the definition of the residuals (3^ (i = 0,1,2). We 
call this choice the reduction procedure 1Z and the exact YFS result J2 n da^ 
is independent of it, if it is done according to the rigorous YFS theory. Our 
choice for TZ is explained in [1]. Finally, we should emphasize that, for efficient 
event generation, it is always desirable to generate a background population of 
events according to a set of distributions da'^ which embody all of the general 
features of Eq. (2.3), but which remove unnecessary details, and to restore the 
exact distributions da^ in (2.3) by rejection methods. In Refs. [1] we follow this 
strategy in constructing YFS2; in addition, several changes of variables are used 
to make the background generation simpler and more efficient from the standpoint 
of CPU time. In this way we have realized the YFS theory for e+e~ — > // + wy 
with /3 , and f3 2l where we should emphasize that /3 2 has only been included 
in the second-order leading-log approximation [6]. 

In the next section, we discuss how we extend YFS2 to more general incoming 
// and // initial and final states, as well as the modifications needed to make the 
program applicable in the SSC energy regime. 

3. YFS2 at SSC energies 

In this section we describe how one extends the YFS2 Monte Carlo program in 
Ref. [1] to realize particle interaction at SSC energies. Such an extension involves 
the introduction of new physics (mainly through a modification of the Born cross 
section), numerical problems (due to the very high energies involved, care is needed 
for the accuracy of the formulas), and certain technical problems associated with 
the Monte Carlo weight rejection method. 
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We begin by discussing the modifications made to introduce the new physics 
at SSC energies. The new program, SSC-YFS2, computes the cross section for the 
interaction 



/(Pl) + 7(P2) 



f'(qi) + Y'(q2)+l(h) + ... + 7(kn) , 



(3.1) 



where / is any lepton or quark. This is still not the most general form of SSC 
interactions, because the incoming fermions are of the same type (identical or 
particle- antiparticle). We are currently generalizing the program to include all 
interactions, and shall report on the results shortly [2] . 

The mass parameters m q used for the quarks are the Lagrangian quark masses. 
We have in mind that the overall momentum transfer in the interactions will be 
large compared to the typical momenta inside the proton. In fact, these quark 
mass parameters should strictly speaking be running masses m q (n), where \x is the 
scale at which they are being probed. Such a running mass effect is well-known and 
is readily incorporated in the program, as the accuracy one is interested in may 
dictate. Thus, with this understanding, further explicit reference to the running 
mass effect is suppressed. 

The interactions realized by YFS2 (Eq. (2.1)) involve only an exchange of 7 and 
Z° in the s-channel. For SSC-YFS2, realizing the more general interaction (3.1), 7, 
Z° and W exchange in the t- and -u-channels, accordingly, had to be introduced. 
This was done by generalizing the Born cross section to include the additional 
channels. Moreover, in the case of quark interactions, a gluon exchange was added 
in all three channels. The running strong coupling constant 
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(3.2) 



(33-2n/)ln{^/(A^ s ) 2 } 



was used, where nj is the number of quark flavors below the energy level [i. In 
our case, rif = 6, and therefore the QCD parameter A.^ s is used. It can easily be 



\ At present, the program cannot handle third-generation fermions. 
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related to the experimentally measured parameter A¥ s = 238 MeV: 



a ms = a ms | Af! j j A f = Af s ( ^- ) . (3.3) 




The masses of the top and bottom quarks were set to = 5 GeV and m t = 
250 GeV, respectively, but the results are little affected by their precise values. 

Certain numerical problems arise at very high energies, because of the very 
small value of all ratios m/y/s, where m is the mass of any interacting particle, 
and \fs = 40 TeV is the energy of the incoming fermions in their center-of-mass 
frame. Certain formulas had to be rewritten so that such small numbers would not 
be ignored by the computer when they should not be; if one is not careful, ratios 
of the form 0/0 appear at various places. Working at SSC energies, however, has 
the advantage that all terms of order m 2 / s or higher can be dropped. The error is 
negligible and leads to a considerable simplification of formulas, and consequently 
to a reduction in computer time. 

Next we discuss the event-generation procedure. To perform the integral for 
the total cross section, we first simplify the form of the differential cross section. 
Thus the exact cross section da is replaced by da', so that the integral J da' can 
be performed analytically. The exact cross section is then computed by rejecting 
events according to their weights, 



" = 5? ■ (3 ' 4) 



Apart from simplicity, we require that da' lead to an efficient generation of events. 
In YFS2, da' was chosen to be a constant. In the present case, this is no longer 
possible, because of the presence of the t-channel. The cross section has a singu- 
larity at t = (pi — q\) 2 = of the form 1/t 2 . To account for the singularity, an 
angle cutoff 9q = 100 mrad is introduced. This is in accord with current detector 



7 



capabilities, and can be changed at will. A crude cross section da' of the form 

da' = A + ^ , (3.5) 

was chosen, where the constants A and B depend on the interaction?" When a u- 
channel also contributes, a similar term of the form 1/u 2 must be added to account 
for the singularity at u = (pi — q^) 2 = 0. 

Finally, we comment on the choice of the reduction procedure which is needed 
for the definition of the arguments of the YFS residuals P i {i = 0, 1, 2), as explained 
in Section 2. The reduction procedure is more delicate in the presence of the t- 
channel, due to the singularity at t — 0. One has to make sure that the weights (3.4) 
do not become uncontrollably large. This is managed by making t as large as 
possible after the reduction. It is not always possible to increase the reduced t 
so that the weight (3.4) remains below the maximum weight. The object of this 
exercise is to minimize the error originating from the tail of the distribution of 
weights above the maximum weight (which is set to 3, but can be changed if so 
desired). This is accomplished by a somewhat involved reduction procedure, which 
is an adaptation of the similar procedure in BHLUMIl.xx [7]. 

This concludes our discussion of the modifications in the YFS2 program nec- 
essary in order to realize interactions at SSC energies. Next, we present some of 
our results. 



X A fictitious photon mass cutoff was also tested, but it turned out to lead to a large weight 
rejection rate. 
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4. Multiple-photon effects at SSC energies 



In this Section we present some results on the effects of multiple-photon initial- 
state radiation on the incoming qq and qq "beams" at SSC energies using our YFS 
Monte Carlo event generator SSC-YFS2 Fortran. Our objective is to determine 
the size of these effects with an eye toward their incorporation into SSC physics 
event generators. This latter step will be taken up elsewhere [2]. 

We consider, as illustrated in Fig. 2 (the kinematics is summarized in the 
figure), 

q(Pl) + ( q\P2) — > q'(qi) + i q ) '(q2) + l(h) + ... + -y(k n ) , (4.1) 

at y^s = 40 TeV for q, q' = u, d, s. For definiteness, we will illustrate our results 
with q = u,d, where we use m u = 5.1 x 1CT 6 TeV, 1714 = 8.9 x 1CT 6 TeV, and view 
y/s = 40 TeV as our worst-case scenario. The more typical [8] value yfs g40 TeV 
6.7 TeV is also presented here for completeness. For these respective input 
scenarios, we shall discuss the following distributions: the number of photons per 
event, the value of v — (s — s')/s, where s' = (qi + q2) 2 , and the squared transverse 
momentum of the outgoing 717 state. These distributions give us a view of the 
effect of this multiple-photon radiation on the incoming quarks and (anti)quarks 
in the SSC environment, where one is really interested in pp — ■> H + X, where H 
is the Standard Model Higgs particle. 

Considering first the number of photons per event, we have the results in Fig. 3. 
There, we show that for the uu incoming beams, the mean number (n 7 ) of radiated 
photons is 0.85 ±0.92 (it is similar for y/s = 6.7 TeV). This should be compared to 
the dd incoming state, where (n 7 ) is 0.21 ± 0.45. For reference, we recall [1] that 
at LEP/SLC energies, the corresponding value of (n 7 ) is, for the incoming e + e _ 
state, ~ 1.5 ± 1.0. Hence, we see here one immediate effect of the high energy of 
the SSC incoming beams: the initial uu-type state will radiate a significant number 
of real photons, with a consequent change in the observed final-state character. In 
particular, the issue of how much energy is lost to photon radiation is of immediate 
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interest, for this energy is unavailable for Higgs production by uu (or dd) and, 
further, it may fake a signal of H — > 77 if we are unlucky. Accordingly, we now 
look at the predicted distribution of 



where s' = (q\ + q^) 2 is the squared invariant final fermion pair mass. If only one 
photon is radiated, v is just the energy of this photon in the center-of-mass system 
of the incoming beams (in units of the incoming beam energy). 

What we find for v is shown in Fig. 4 for the uu — > uu + ny case (the dd — ► 
dd + n7 case is similar). We see the expected shape of v from Ref. [1], and its 
average value is (v) = 0.05 ± 0.09. Hence, ~ 10% of the incoming energy is 
radiated into photons; this energy is not available for Higgs production and hence 
it is crucial to fold our radiation into the currently available SSC Higgs production 
Monte Carlo event generators [9] and to complete the development of our own YFS 
multiple-photon (-gluon) Higgs production Monte Carlo event generator, which is 
under development and will appear elsewhere [2]. 

Given that we know we have, in the SSC environment, significant multiple- 
photon radiation effects, the question of immediate interest is how often the trans- 
verse momenta of two photons are large enough that they could fake a H — > 77 
signal. We will answer this very important question in detail in the not-too-distant 
future when our complete Higgs production YFS Monte Carlo event generators are 
available [2]. However, here we can begin to study this question by looking into the 
transverse momentum distribution of our YFS multiple-photon radiation in, e.g., 
uu — > uu + ny. This is shown in Fig. 5, where we plot the total transverse momen- 
tum distribution of the respective YFS multiple-photon radiation. What we find 
is that , for y/s = 40 TeV, the average value of this total transverse momentum is 
(in the incoming uu center-of-mass system) 




(4.2) 




(4.3) 
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where k{ {% — 1, . . . , n) are the four- momenta of the n photons. (For ^/s = 6.7 TeV, 
this average is (0.0186 ± 0.0136) v /i.) Hence, for the SDC acceptance cut of 
^| lntan(#/2)| = \t)\ < 2.8, or 6{ > 122 mrad, this means that there may be 
some possible background to H — > 77 for, e.g., mjj m 150 GeV. Such effects will 
be discussed in detail elsewhere [2] . 

Finally, concerning the overall normalization, we find that the Born cross sec- 
tion is corrected according to the results in Table 1. This shows clearly that the 
higher-order effects change the normalization by approximately 5%. This sets the 
level at which precise simulations of SSC physics must take the higher-order effects 
from multiple photons into account. 

5. Conclusions 

We conclude that our initial study of YFS multiple-photon radiation in the SSC 
physics environment shows that any Monte Carlo event generator which hopes to 
achieve an accuracy of order 10% in the SSC physics simulations must treat the 
respective effects in a complete way. In this paper, we have computed these effects 
for incoming quark- (anti) quark states at SSC energies using the Monte Carlo event 
generator SSC-YFS2 Fortran based on our original YFS2 Monte Carlo in Ref. [1]. 

Specifically, using our SSC-YFS2 Monte Carlo event generator for q { q ) — > 
+ Try, at y/s = 40 TeV, we find that for an initial uu state, the mean num- 
ber of radiated photons is 0.85 ± 0.92, so that the multiple-photon character of 
the events must be taken into account in detailed detector simulation and physics 
analysis studies. Further, the mean value of v — (s — s')/s is 0.05 ± 0.09 and the 
average total squared transverse momentum (k\ to ^j is 0.025 ±0.002 s. Hence, the 
impact of these event characteristics on Higgs production in general and on the 
H — > 77 scenario in particular must be assessed in detail. Such assessment will 
appear elsewhere [2]. 

In conclusion, we can say that the initial platform for precision SSC elec- 
troweak physics simulations on an event-by-event basis using our YFS Monte Carlo 



11 



approach [1] has now been established. We look forward with excitement to its 
complete development for all such electroweak phenomena and to its extension to 
the SSC QCD processes as well [2]. 
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FIGURE CAPTIONS 



1) The process e+(pi) + e (p 2 ) -> /(<7l) + /(<?2) + 7(^1) + • • • + l{K)- 

2) The SSC process q( Pl ) + i q ) (p 2 ) -> g'(gi) + <g>'(g 2 ) + 7(^1) + • • • + 7(*n), 
where g, g' = u, d, s. 

3) Histogram of the photon multiplicity in uu — > mm + 717 for |r/| < 2.8: 

(aj v^i = 40 TeV; (b) y/s = 6.7 TeV. Here, v m i n = 10~ 6 -we have shown in 
Ref. [1] that the cross section does not depend on v m i n . 

4) ^-distribution for uu —> uu + 717, where v — (s — s')/s and s' = (q\ + g 2 ) 2 is 
the squared final uu invariant mass. Here, \rj\ < 2.8: (a) \/s = 40 TeV; (b) 
y/s = 6.7 TeV. 

5) Total transverse momentum distribution of the photons in uu — » uu + wy for 
\r]\ < 2.8 in units of s: (a) ^1 = 40 TeV; f&j = 6.7 TeV. 

TABLE CAPTIONS 

1: Sample output for uu — > mt + 717 at ^/s = 40 TeV and |r/| < 2.8. The entries 
in the table are largely explained therein: XSEC = cross section, WT = 
event weight, and BORN = Born cross section. 
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